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We parameterize the recently proposed generalized Skyrme effective force (GSEF) containing 



^ . 
.^ ■ extended density dependence. The parameters of the GSEF are determined by the fit to several 

^^ , properties of the normal and isospin-rich nuclei. We also include in our fit a realistic equation of 

state for the pure neutron matter up to high densities so that the resulting Skyrme parameters 

iy~, ■ can be suitably used to model the neutron star with the "canonical" mass (~ 1.4Mq). For the 

(•~^ . appropriate comparison we generate a parameter set for the standard Skyrme effective force (SSEF) 

CD using exactly the same set of the data as employed to determine the parameters of the GSEF. We 

^D ' find that the GSEF yields larger values for the neutron skin thickness which are closer to the 

recent predictions based on the isospin diffusion data. The Skyrme parameters so obtained are 

O I employed to compute the strength function for the isoscalar giant monopole, dipole and quadrupole 



resonances. It is found that in the case of GSEF, due to the the larger value of the nucleon effective 
mass the values of centroid energies for the isoscalar giant resonances are in better agreement with 
the corresponding experimental data in comparison to those obtained using the SSEF. We also 
present results for some of the key properties associated with the neutron star of "canonical" mass 
and for the one with the maximum mass. 

PACS numbers: 21.10.-k,21.65+f,24.30.Cz,21.60jz,26.60.+c 



I. INTRODUCTION 

The density dependent Skyrme type effective nucleon-nucleon interaction jjj within the 
Hartree-Fock (HF) approximation has been one of the most successful and popular micro- 
scopic tools to describe the ground state properties of the finite nuclei as well as that of 
the symmetric nuclear and pure neutron matters. The pioneering work ^ of implementing 
the Skyrme type effective force having only a linear density dependence was carried out to 
reproduce the experimental data on the binding energy and charge rms radii. The linear 
density dependence gave rise to the value of nuclear matter incompressibility coefficient K^o 
in the range of 300 — 400 MeV which is much higher compared to the experimental value 
~ 220 MeV. It was proposed in Ref. J3| that in order to obtain a reasonable value of K^ 
the linear density dependence in the Skyrme effective force must be modified to p" with a 
lying in between | to |. Since then, numerous attempts have been made to parameterize 



the standard Skyrme effective 



brce (SSEF) which contains only a single density dependent 



term (e.g., see Ref. |^, y, M, LD, Isl, Isl, llil lui ll^, lld|)- Recently several parameter sets for 
the SSEF are obtained by fits to large set of data comprising the binding energy, charge 
rms radii and single-particle energies for the nuclei ranging from normal to isospin-rich ones. 
Further, the Skyrme parameters are constrained by demanding that the equation of state 
(EOS) for the pure neutron matter should be reasonable up to the densities relevant for 
studying the properties of neutron star J9|, |l3, ll3| . 

Recently, studies involving the generalized Skyrme effective force (GSEF) have been re- 
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visited jU, |15|, |16|, |17(. Generalization of the Skyrme effective force can be realized by 
adding several density dependent terms to each of the three, namely, local, non-local and 



m 



are 



spin-orbit parts of the SSEF. The multi density-dependent terms used in Refs 
of the form p'^^^ with v = 1,2, and 3 connected to the Fermi momentum [kp] expansion 
of the Brueckner G matrix. Similar density dependence is also found in the energy density 
functional derived within the chiral perturbation theorv containing the contributions from 
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the GSEF con- 



one and two pion exchange diagrams up to three- loops |lq]. In Refs. 
taining extended density dependence in the local part of the Skyrme interaction is used only 
to study the properties of nuclear matter. It is found that the extended density dependence 
gives rise to reasonable EOS for the infinite nuclear matter for any given asymmetry. It has 
been also demonstrated that unlike in the case of SSEF the value of K^o and the isoscalar 



nucleoli effective mass m* can be determined independently using the GSEF. Thus, one 
may expect the GSEF to reproduce concurrently the experimental values for the i^oo and 
fi^article energies in nuclei. We would like to mention that some early attempts 
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22i | to generalize the density dependence of Skyrme t ype effective interaction 
were significantly different compared to the one proposed recently [15|,|l6|. In Ref. j22j the 
exponents of the density dependent terms are obtained by a fit to several nuclear observ- 
ables. Moreover, most of these forces are plagued by an undesirable feature that they give 
collapse in the nuclear matter at high densities which leave them unsuitable for their use for 
the study of neutron star. 

In the present paper we determine the parameters of the GSEF by the fit to a set of 
experimental data for the binding energies, charge rms radii, single particle energies and 
rms radii of valence neutron orbits. Our data set used in the fit consists of 13 spherical 
nuclei, namely, ^^O, ^^O, ^°Ca, ^^Ca, ^^Ni, ^gni, esNi, ^s^i^ ssg^^ 902^, loogj^^ i32gj^ ^^^^ 

^°^Pb. We also fit the experimental data for the breathing mode energies for the ^°Zr, and 
208p|-| jiuclei. Further, we constrain the value of Skyrme parameters by including in the 
fit a realistic EOS for the pure neutron matter. The chi-square minimization, required to 
obtain the best fit parameters, is achieved by using the simulated annealing method (SAM) 
as recently implemented jl3| to determine the parameters of the SSEF. 

The present paper is organized as follows. In Sec. |n]we briefly outline the form of the 
GSEF and the corresponding energy density functional adopted in the present work. In this 
section, we also mention in short the strategies used to evaluate center of mass corrections 
to the binding energy and charge radii. In Sec. II I II we briefly describe a procedure for mini- 
mization of the x^ function based on the SAM and present the set of the experimental data 
along with the constraints used in the fit to determine the values of the Skyrme parameters. 
In the same section we list the values of the parameter sets for the GSEF and SSEF. In Sec. 
II VI we present our results for the three different fits carried out in this work. We also present 
results for the isoscalar giant monopole, dipole and quadrupole resonances and some key 
properties of the neutron stars obtained using the newly generated parameter sets. Finally, 
in Sec. Mwe summarize our main results. 



II. GENERALIZED SKYRME EFFECTIVE FORCE 



The GSEF used in Refs. Il5l Il6|| can be written as 
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with z = 1, 2, 3, ... and Oj = i/3. In Eq. (P), r = ri — r2, P 



2i ' 



(1) 

P' is complex conjugate 
of P acting on the left and a = ai -\- a2, P^ = | (1 + ai • (T2). The SSEF can be obtained 
from Eq. (^ simply by setting t^i = for i ^ I and ai is normally taken^o be less than 
unity. For instance, in the SIII force 



23| «! is unity and for the SLy forces J9| ai = 1/6 was 
used. It may be noted that we have considered the extended density dependence only for 
the local term in Eq. (0). However, on the same analogy the non-local and spin-orbit terms 
can also be extended to have density dependence. 
The total energy E of the system is given by. 



E = / nir)d'r 



(2) 



where, 7i(r) is the Skyrme energy densi 
the time-reversal invariance is given by 



'unctional corresponding to Eq. ((H) which under 
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JC + Hs + 'Hp + HcS + Tifin + Tiso + '^sg + '^Coul 
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where, /C = 2^t is the kinetic energy term, Hs is the zero-range term. Tip the density 
dependent term, Hes an effective-mass term, Tifin a finite-range term, Tiso a spin-orbit term, 
Tisg a term due to tensor coupling with spin and gradient and T^coui is the contribution to 
the energy density for protons due to the Coulomb interaction. For the Skyrme interaction 
of Eq. ((T]), we have, 



Hs = -to [(2 + xo)p' - (2x0 + l)(pj + p; 



(4) 



^P = i E^3.P"' [(2 + X3.)p' - (2x3. + l)(pj + pD] 



(5) 
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Kff = ^ [ti(2 + xi) + h{2 + X2)] TP + - [t2{2x2 + 1) - ti(2a;i + 1)] (r^pp + r„p„), (6) 



-^ [3ti(2xi + 1) + t2(2x2 + 1)] [(Vpp)2 + (Vp„)2] , (7) 

Wo 

^so = ^ [J ■ Vp + Jp • Vpp + Jn ■ Vpn] , (8) 

Hsg = -J^l^ia^l + ^2X2)J' + Y^(tl - t2) [Jp' + Jn'] , (9) 

ncAr) - \e^pAr) j '-^ - \e^ P^r) (^) ^'^ ■ (10) 

Here, p = Pp + pn, t = Tp + Tn, and J = Jp + J„ are the particle number density, kinetic 
energy density and spin density with p and n denoting the protons and neutrons, respectively. 
We have used the value of h'^/2'm = 20.734 MeV fm^ in our calculations. We would like to 
emphasize that we have included the contributions from the spin-density term as given by 
Eq. Q so that the corresponding contributions can be considered consistently in evaluating 
the Landau parameter G'q ^j. Although the contributions from Eq. Q to the binding 
energy and charge radii are not very significant, but they may contribute significantly to the 
value of the Landau parameter G'q. 

The single-particle wave functions (pi and corresponding single-particle energies Ei are 
obtained by solving the HP equations given by. 
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(f)i{r,q) =ei(j)i{r,q), (11) 



where, m* is the effective nucleon mass, Ug and Wg are the central and spin-orbit parts of 
the mean field potentials. The expressions for the m* and Wg can be found in Ref . p| . The 



f/,(r) = -to[(2 + a;o)p(r)-(l + 2a;o)p,(r)] 



expression for Ug is given as 
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-4 [3ti(2 + xi) - t2(2 + 0:2)] VV(r) 
lo 

+ 4 [3ti(l + 2xi) + t2(l + 2x2)] V%(r) 
lb 

-^W^o (V ■ J + V ■ J,) + 5,,p\/coui(r), (12) 

where, q = p,n. 

To this end we briefly outhne the corrections made to the HF resuhs for the binding 
energies and charge rms radii arising from various effects as follows. The center of mass (CM) 
corrections to the binding energies and charge rms radii are made using simple prescriptions 
discussed in Ref . Il3l UB . We also consider the corrections to the charge rms radii due to the 

nn 

spin-orbit effect as well as the charge distributions of the neutron and proton J9|, |26|. The 
contribution from the Coulomb exchange term in Eq. |T0|) is dropped in order to compensate 
for the effects of the long-range correlations J22| on the Coulomb displacement energy (CDE) 
in mirror nuclei. We also add the modified Wigner term 

N-ZV] ,, ,,. „, f /A^^" 



E^ = V^exp{-Xl-^—j \+V^,\N-Z\exp^-{ — \ } (13) 

to the binding energy. Values of the parameters \4), V^', X and Aq are determined by the 
fit. 



III. PARAMETERIZATION OF THE GENERALIZED SKYRME EFFECTIVE 
FORCE 

The parameters of the GSEF can be determined by minimizing the value of x^ which is 
given as, 






^ =^r^wT.[ .. (14) 



where, Nd is the number of experimental data points and Np the number of fitted param- 
eters. The Oi stands for theoretical error and M^^^ and Mf^ are the experimental and 



the corresponding theoretical values, respectively, for a given observable. The values of x^ 
depends on the Skyrme parameters, since, the M*^ in Eq. (fT^ is calculated using the HF 
approach with a Skyrme type effective nucleon-nucleon interaction. Thus, it is clear that the 
best fit parameters can be obtained using an appropriate method for the x^ minimization 
together with a set of experimental data. We briefly describe below the SAM for the x^ 
minimization and the set of experimental data selected to obtain the best fit parameters. 

A. SAM algorithm for x^ minimization 

The concept of SAM is based on the manner in which liquids freeze or metals recrystallize 
in the process of annealing. In an annealing process a metal, initially at high temperature and 
disordered, slowly cools so that the system at any time is in a thermodynamic equilibrium. 
As cooling proceeds, the system becomes more ordered and approaches a frozen ground 
state at zero temperature. The SAM is an elegant technique for optimization problems of 
large scale, in particular, where a desired global extremum is hidden among many local 
extrema. This method has been found to be an extremely useful tool for a wide variety of 
minimization problems of large non- linear systems in many different areas of science (e.g.. 



see Refs 



a 



30 



3l|). Very recently 



32, |33|, the SAM was used to generate some initial 



trial parameter sets for the point coupling variant of the relativistic mean field model. 

In Ref. lJ| we have implemented the SAM to the problem of searching the global 
minimum in the hypersurface of x^ function as defined by Eq. ()14|). Here too we shall 
use the SAM to determine the parameters of the GSEF. In the SAM one needs to specify 
the appropriate annealing schedule together with the parameter space (i.e., the range of the 
values of the parameters) in which the best fit parameters are to be searched. Similar to that 

n 

in Ref. jl3|, in the present work we have employed a moderately faster annealing schedule 
such as Cauchy annealing schedule given by, 

T{k) = T,/k (15) 

where, Tj is the initial value of the control parameter (also viewed as effective temperature) 
and T{k) with k = 1,2, 3, ... is the control parameter at k—th step. The value of k is increased 
by unity after lOOA'p reconfigurations or lOA'p successful reconfigurations whichever occurs 
first. The value of Tj is taken to be 1.25 which is same as that used in Ref. [l3|. We keep 



on reducing the value of the control parameter using Eq. (P3|) in the subsequent steps until 
the effort to reduce the value of x^ further becomes sufficiently discouraging. In Ref. [l3| 
instead of the range of the values of the Skyrme parameters we used range of the values 
of the quantities associated with the symmetric nuclear matter. Since, it was possible to 
identify the number of nuclear matter quantities linearly related to the equal number of 
parameters for the SSEF. But, in the case of GSEF the numbers of parameters are larger 
and we shall define the parameter space directly in terms of the range of the values each 
of the Skyrme parameter can take. In Table H] we give the lower and upper limits for the 
values of the Skyrme parameters denoted by vq and vi, respectively. The quantity d in the 
penultimate column denotes the maximum displacement allowed in a single step for a given 
Skyrme parameter during the reconfiguration process (see also Ref. [l3| for detail). The last 
column labeled as Vj„ contains initial values for the Skyrme parameters used as a starting 
point for the SAM. These values for the parameters are the same as that for SGI force 
^ which has p^'^ density dependence. The parameters of the Wigner term (Eq. IT3|l are 
taken from Ref. J2q]. The lower and upper limits and the initial values for the parameters 
ai are taken to be the same indicating that these parameters are kept fixed during the x^ 
minimization. 

B. Data used in the fitting procedure 

We now summarize our selection of data and corresponding theoretical errors adopted 



in the 



itting procedure. In Ta ble ITTl we list our choice of the data along with their sources 
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. It must be noted that in addition to the data on the 
binding energy, charge radii and single particle energies, the values of the Skyrme parameters 
are further constrained by including in the fit the experimental data for the radii of valence 
neutron orbits, breathing mode energies together with a realistic EOS for the pure neutron 
matter upto the densities (~ 0.8 fm~^) relevant for the study of neutron stars. For the 
binding energy we use in our fit the error of 1.0 MeV except for the ^°°Sn nucleus. The 
binding energy for the ^°'^Sn nucleus is determined from systematics and are expected to 
have large errors. Thus, we assign it a theoretical error of 2.0 MeV. For the charge rms radii 
we use the theoretical error of 0.02 fm except for the case of ^^Ni nucleus. The charge rms 
radius for the ^^Ni nucleus is obtained from systematics and we use the theoretical error 



of 0.04 fm. For the rms radii of the valence neutron orbits in ^^O and ^^Ca nuclei we use 

371 138| respectively. The theoretical error 



rt,(z/l(i5/2) = 3.36 fm and r.a{i'lfT/2) = 3.99 fm, 
taken for the rms radii for the valence neutron orbits is 0.06 fm. We must point out that the 
choice of the theoretical error on the rms radii for the valence neutron orbits is due to the 
large uncertainties associated with their extraction from the experimental measurements. To 
be consistent with the way these valence neutron radii are determined, we do not include the 
center of mass correction to these data. For each of the 22 single-particle (S-P) energies in 
the 208p]^ nucleus, we have used the theoretical error of 1.0 MeV in our fit. The experimental 
data for the breathing mode constrained energies En included in our fit are 17.81, and 14.18 
MeV fo. the «'Zr and -Pb nude, Q. respective... with the theo.et.ea. e^o. taken to be 



SO include in the fit 
We use 15 data 



0.5 MeV for the ^°Zr nucleus and 0.3 MeV for the ^o^Pb nucleus. We a 

the EOS for the pure neutron matter of a realistic UV14+UVII model 

points for the EOS corresponding to the densities up to 5po- The theoretical error on each 

of these 15 data points are taken to be 2.0 MeV. 

C. Parameters for the generalized Skyrme effective force 

Following the fitting procedure described in Sec. IIII Al together with the list of data given 
in Table ITTl we have generated two different parameter sets for the GSEF. These parameter 
sets are given in Table IIIII The parameter set GSkI and GSkll are obtained as follows; 
GSkI includes the density dependent terms proportional to p"' with ctj = 1/3, 2/3 and 1, 
whereas, GSkll force has only p^^^ and p^/^ density dependent terms. The exponents for 
the density dependence are taken to be 1/3 and 2/3 for the GSkll interaction due to the 
following reasons. The p^^^ density dependence decouples the effective mass and the nuclear 
matter incompressibility coefficient 1^. On the otherhand, the SSEF favours the exponent 
of density dependence to lie in between 1/3 — 1/6 in order to give an acceptable value for 



the nuclear matter incompressibility coefficient. In the last column of Table UTTl we give the 
values of parameters corresponding to the SSEF 9]. We name this force as SSk which is 
obtained using the same set of data as those for the GSkI and GSkll forces. We shall see in 
the following section that the quality of the fit for the GSkI force is better in comparison to 
the ones obtained using GSkll and SSk forces. It can be seen from Eq. (Q) that the density 
dependence in the GSkI force requires six parameters, namely, tsi and xsi with i = 1,2 and 



3, whereas, in the SSk force the density dependence is specified by three parameters, t^i, X31 
and ai. Thus, fit to the GSEF considered here requires three additional parameters. 

IV. RESULTS AND DISCUSSIONS 

We generate the parameter sets GSkI and GSkll for the GSEF using the experimental 
data on the bulk properties for the nuclei ranging from the normal to isospin-rich ones. The 
values for the Skyrme parameters are further constrained by including in fit the experimental 
data on the breathing mode energies and a realistic EOS for the pure neutron matter. For 
the appropriate comparison we also obtain a parameter set named SSk for the SSEF using 
the same set of data as used to determine the parameters of the GSEF. The values of 
parameters for all the three forces in consideration are given in Table IIIII In the following 
subsections we shall present the results obtained using these three different parameter sets. 

A. Infinite symmetric and asymmetric nuclear matter 

In Table ITVI we give the values of various quantities associated with the symmetric nuclear 
matter at the saturation density obtained for GSkI, GSkll and SSk forces. We see that the 
values for all the quantities except for the m* and 

given in Table IIVI are quite close to each other for all the three forces. The values of m* 
and L for the GSkI and GSkll forces are little larger than that for the SSk force. The value 
of m* ~ 0.8m for the GSkI and GSkll forces is highly desirable in order to appropriately 



reproduce the location of the isoscalar giant quadrupole resonance 43|. Further, we shall 
see below that larger value of m* gives rise to the single-particle energies which are in better 
agreement with the experimental data. The values of L at the saturation density obtained 
for the GSkI and GSkll forces in comparison to the one for the SSk force are in better 
agreement with the recently extracted value of L = 88 ± 25 MeV J4^ . 

We consider the behavior of the symmetr y energ y coefficient J{p) for densities relevant to 
the study of neutron stars. It is well known (45l. l46l| that the values of J{p) and the resulting 
EOS for the pure neutron matter at higher densities (p > 2po, po = 0.16 fm~^) are crucial 
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in understanding the various properties of neutron star. For example, the proton fraction at 
any density depends strongly on the value of J{p) at that density, which in turn affects the 
chemical compositions as well as the cooling mechanism of the neutron star |47|. Yet, no 
consensus is reached for the density dependence of J{p)- We display in Fig. d our results 
for the variation of the symmetry energy J as a function of the nuclear matter density p. 
We see that the value of J increases with density. The slope of J{p) being positive upto 
3po for GSkI, GSkll and SSk forces indicate that these interactions can be used to study 
the properties of the neutron star with masses around IAMq J48|. In Fig. El we plot the 
EOS for the nuclear matter for various proton fraction ( Yp = Pp/p) as a function of density. 
The solid squares represent the EOS for the pure neutron matter {i.e.,Yp = 0) for a realistic 
UV14+UVII model j423. It can be seen that the EOS for the pure neutron matter obtained 
using GSkI force is in better agreement with the one obtained for the UV14+UVII model. 

B. Fit to the nuclear bulk properties 

In Fig. El we present our results obtained using the GSkI, GSkll and SSk forces for the 
relatives errors in the values of total binding energies for the nuclei used in the fit. The rms 
error for the total binding energy are 1.18, 1.72 and 1.28 MeV for the GSkI, GSkll and SSk 
forces, respectively. Thus, We can say that the quality of the fit to the binding energies 
for the GSkI and SSk forces are more or less the same and they are much better than that 
for GSkll force. The binding energy difference i?(^^Ca) — i?(^^Ni) = 66.82, 66.57 and 66.95 
MeV for the GSkI, GSkll and SSk interactions, respectively, compared to the experimental 
value of 68.85 MeV. The said difference for the SKX interaction |lO| is 66.3 MeV which 
is about 0.5 MeV lower as compared to our present results. On the other hand, most of 
the Skyrme interactions which include the contribution from the exchange Coulomb term 
yield -B(^^Ca) — i?(^^Ni) ^ 63 MeV, which is about 6 MeV lower than the corresponding 
experimental value. To this end we would like to remark that the Wigner corrections (Eq. 
(I13|) ) are more significant for the GSk forces as compared to the one for the SSk force. It 
may be verified by using the values of the Wigner parameters given in Table UTTl that for the 
symmetric nuclei the Wigner correction is greater than 2.5 MeV for the GSk forces which 
is more than twice than that for the SSk force. Also for asymmetric nuclei we find that the 
GSk forces has larger Wigner corrections for the ligher nuclei. For example, in case of ^^O 
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the Wigner correction is about 4.0(0.5) MeV for the GSk(SSk) forces. In Fig. I^we present 
our results for the relative errors in the values of charge rms radii for the nuclei used in the 
fit. The rms error for the charge rms radii are 0.023, 0.026 and 0.021 fm for the GSkI, GSkll 
and SSk forces, respectively. The quality of fit for the charge rms radii is more or less same 
for all the three forces considered here. 

In Fig. El we give the values for the single-particle energies for the ^°^Pb nucleus used 
in the fits. We see that the values of the single-particle energies for the case of GSkI and 
GSkll forces are in slightly better agreement with the experimental data as compared to 
the ones obtained for the SSk force. This is mainly due to the weakening of the correlations 
between the effective mass m* and the incompressibility coefficient K^ in the case of the 
GSEF [l5|, |l6| . In Table \V\ we present our results for the rms radii for the valence neutron 
orbits and the breathing mode constrained energies used in the fit. The breathing mode 
constrained energy Eq is given by. 



E. = J^ (17) 

where, fully self-consistent values for the energy- weighted (mi) and inverse energy- weighted 
(?7i_i) moments of the strength function calculated using random phase approximation 
(RPA) for the breathing mode (or isoscalar giant monopole resonance) can be obtained 
using the double commutator sum rule and the constrained HF method |l3|, |49|, respec- 
tively. We see that all the Skyrme forces under consideration gives the similar fit to the 
experimental data presented in Table |3 

C. Neutron skin thickness 

The neutron skin thickness r„ — Vp is the difference between the rms radii obtained 
using the density distributions for the point neutrons and point protons. It is well known 
pOL l5l| that the accurate measurement of the neutron skin thickness would place a stringent 
constraint on the density dependence of the nuclear symmetry energy. In particular, it has 
been shown [44] that the neutron skin thickness in heavy nuclei are linearly correlated with 
the slope of the symmetry energy coefficient L. The value of L is poorly known due to 
the large uncertainty associated with the measurement of the neutron skin thickness. The 
measured values of the neutron skin thickness is very much probe dependent. Current data 
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indicate tha t the neutron skin thickness in the ^'J^Pb nucleus he in the range of 0.10 
0.28 fm 



53]. The proposed experiment at the Jefferson Laboratory on parity violating 
electron scattering from ^°'^Pb is expected to give another independent and more accurate 
measurement (wit hin 0.05 fm) of its neutron skin thickness. Recent analysis of the isospin 



diffusion data |4^ yields L = 88 ± 25 MeV which restricts the values of the neutron skin in 
the ^^^Sn and ^o^Pb nuclei to 0.29 ± 0.04 fm and 0.22 ± 0.04 fm, respectively. In Fig. Hwe 
present our results for neutron skin thickness obtained using GSkI, GSkll and SSk forces. 
The values for neutron skin thickness for the GSkI and GSkll forces are larger than those 
for the SSk force. It is because of the fact that the values of L for the GSkI and GSkll forces 
are larger than that for the SSk force (see Table HV|) . The values of neuron skin thickness 
obtained using the GSk forces for the ^^^Sn and ^ospb nuclei are in better agreement with 
the recent predictions 0.29 ± 0.04 and 0.22 ± 0.04, respectively. 

D. Isoscalar giant resonances 

We consider now the results for the isoscalar giant monopole, dipole and quadrupole 
resonances obtained using the GSkI, GSkll and SSk forces. The centroid energ y fo r the 
isoscalar giant monopole resonance (ISGMR) depends mainly on the value of K^o J53|. The 
centroid energy for the overtone mode of the isoscalar giant dipole resonance (ISGDR) is 
governed by the values of the m* as well as K^o [54], whereas, the centroid energy of the 
isoscalar giant qua drup ole resonance (ISGQR) being a surface oscillation depends merely 



on the value of m* 43| . As we mentioned earlier, the correlations between the m* and Kqq 
is weaker for the GSEF which may give rise to a better agreement of the experimental data 
for the isoscalar giant resonances. 

The excitation energy dependence of the strength function S{E) for the isoscalar giant 
resonances is obtained using the RPA Green's function G as J55l | . 

5(i?) = V|(0|F|n)|^5(E-E„) = -Im[Tr(/G/)], (18) 

^ — ^ TT 

n 

where, 

A 

F = Y. />^) (19) 
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is the scattering operator with 

r^Yooir) ISGMR 



/(r) 



r 



3 



rir)Yioir) ISGDR (20) 



r^Y2o{f) ISGQR 

for the various isoscalar giant resonances considered. The quantity t] in Eq. (j2(Jj) for the 



case of ISGDR is taken to be equal to |(r^) in order to ehminate the spurious contributions 
arising from the center of mass motion. The RPA Green's function appearing in the Eq. 
(fTS|) is obtained by discretizing the continuum 5y|. The continuum is discretized by a box 
of 30 fni with the radial mesh size of 0.3 fm. The strength functions are smeared using 
a Lorentzian with the smearing parameter taken to be equal to 1.0 MeV for the purpose 
of plotting. We must mention that the contributions due to the spin-orbit and Coulomb 
interactions are ignored in evaluating the RPA Green's function, but, they are present in 
the HP calculations. However, neglect of spin-orbit and Coulomb interactions do not have 
any important bearing on our findings. It has been shown in Refs. (57L l58l | that ignoring 
the spin-orbit and Coulomb interactions in the RPA calculations do not affect the values 
of the centroid energies for ^°Zr and ^°^Pb nuclei as the effects of ignoring the spin-orbit 
interaction is counterbalanced by that for the Coulomb interaction. Further, in the present 
work we are mainly interested in the differences between the behaviour of the isoscalar giant 
resonances for the different interactions. 

In Figs. EK -Eb we have plotted the strength function for the ISGMR, ISGDR and ISGQR 
for the ^ospb nucleus obtained using GSkI, GSkll and SSk forces. We see that for the case 
of ISGMR all the three forces give similar results. But, for the case of ISGDR and ISGQR 
the peaks at the high energy for the GSkI and GSkll forces are lowered by about 0.5 MeV in 
comparison to those for the SSk force. Similar trends are also observed in the ^^Zt nucleus 
(not shown). In Table IVTl we present our results for the centroid energy ( ^ ), constrained 
energy J^ and scaling energy J^ for the ISGMR, ISGDR and ISGQR for the ^^Zr and 
^°*Pb nuclei obtained using GSkI, GSkll and SSk forces. The values of the k-th moment r/ifc 
are obtained using the strength function (Eq. IT^ as, 

mk= [ ' E^S{E)dE (21) 

where, Ei and E2 are the minimum and maximum values of the excitation energy over which 
the integration is performed. The range of the excitation energy Ei — E2 used to obtain 
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various energies for the isoscalar giant resonances presented in Table IVII are as follows. 
To evaluate the ISGMR energies we have used Ei — ii^2 to be — 40 MeV for all the 
cases considered. For obtaining the values of ISGDR and ISGQR energies for the ^"Zr 
(^°^Pb) nuclei we use the integration range for the excitation energy as 18 — 40(16 — 40) and 
9.5 — 40(7.5 — 40) MeV, respectively. The non-zero values for the Ei for the case of ISGDR 
and ISGQR are so choosen that the low-lyings peaks are excluded from the calculation of the 
k—th moments given by Eq. ^IT\ . The strength functions used in Eq. ()2ip to evaluate the 
moments rrik are smeared using a small value of the smearing parameter r/2 = 0.05 MeV. 
It is clear from this table that the centroid energies for the ISGDR and ISGQR obtained 
using GSkI and GSkll forces are in better agreement with the corresponding experimental 
data in comparison to those obtained using SSk force. 



E. Properties of the neutron star 

In this section we present our results mainly for the non-rotating neutron star obtained 
using the newly generated Skyrme parameter sets. In particular, we have studied the rela- 
tionship between the mass and radius for the non-rotating neutron star and obtained some of 
the key properties for the neutron star with the "canonical" mass of 1.4M0 and for the one 
with maximum mass Mmax- Bulk properties of spherically symmetric non-rotating neutron 
star can be obtained by integrating the Tolman-Oppenheimer-Volkoff (TOV) equation given 
as. 



dP 
dr 



Ge{r)M{r) 



[l 1 ''^^^l 




L ' e[r) \ 



and 



dM{r) 



1 + 


47rr3p(r) 


M[r)c^ 


47rr 


\{r) 



2GM{r) 



rc^ 



dr (? 



(22) 



(23) 



where, the quantities P, e and M are the pressure, energy density and the gravitational 
mass of the neutron star, respectively, which depend on the radial distance r from the 
center. Numerical solutions of the Eqs. (|22j) and (|23p can be easily obtained for given 
initial values for the pressure -P(O) at the center and using M(0) = 0, provided, the relation 
between P and e is known. For lower densities, we use Baym-Pethick-Sutherland EOS 59|, 
matching onto the Skyrme EOS at p ~ 0.5po and going down to p = 6.0 x 10^^^ fm~^. At 
densities larger than 0.5po we use the Skyrme EOS which is obtained by adding the nuclear 
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Pp 


= Pe+ PtM, 


Pn 


= /ip + /ie; 


Pe 


= Pf.- 



and leptonic contributions . The nuclear part of the energy density is obtained using Eq. 
and following the Ref. 60|. In the leptonic sector we have considered the contributions 
from the electrons and niuons with the corresponding energy densities Ee and e^ given by, 

ei = ^ r k^^k^ + m^dk, (24) 

where, / = e (electron) or /i (muon) in the above equation and kK is the corresponding fermi 
momentum. For a given baryon density p, the values of the fermi momenta for the neutrons, 
protons, electrons and muons can be obtained by ensuring that the conditions for charge 
neutrality and chemical equilibrium are satisfied, i.e., 

(25) 
(26) 
(27) 

In Eqs. fEBj) and (fTTj) . fij denotes the chemical potential with j = n,p,e or /i. Thus, the 
energy density depends only on the baryon density and the pressure can be obtained as, 

P = /%^. (28) 

dp 

We now consider our results for the properties of neutron stars for the GSkI and SSk 

forces. Results for the neutron star with the canonical mass of IAMq obtained using the 

GSkI and GSkll forces are very much similar. Because, for both of these forces the density 

dependence of the symmetry energy and the pure neutron matter are quite close up to the 

densities 3po- However, the GSkll force is not appropriate for the study of the neutron stars 

with masses larger than the canonical ones due to the undesirable feature for the symmetry 

energy at densities higher than 3po (see Fig. [T)). In Fig. |Hlwe have plotted the variations 

of the proton fraction {Yp = Pp/p) and electron fraction [Yg = Pe/p) as a function of baryon 

density. The values for Yp and Y^ or Y^ at all the densities considered are not large enough 

to allow the direct Urea process. In Fig. El we display our results for the variations of the 

neutron star mass as a function of the central density pc for the baryons. The relationship 

between the mass of the neutron star and its radius is shown in Fig. ^1 In Table IVllI 

we present our results for some of the key properties associated with the neutron star with 

masses IAMq and Mmax- These properties are quite similar to the ones obtained using 

SLy4 interaction 48|. The quantity Ec/c^ is the energy density at the center. The values 
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for baryon number (A) and binding energy (i?bind) for neutron star given in Table IVIII are 
obtained as follows, 

Jo (l_2GM^l/2• 
and 

^bind = {Am^ - M)^ (30) 

where, rriQ = 1.66 x 10~^^g is the mass per baryon in ^^Fe. We also give in Table IVTTl the 
value of gravitational redshift of photons (^surf) emitted from the neutron star surface which 
is given as. 

In Eqs. (j2ni) - (j^T|) - R is the radius of neutron star with the gravitational mass M. We 
must point out that softening of the EOS at higher densities due to appearance of the 
hyperons and kaons might alter the results. In particular, results for the neutron star with 
the maximum mass is expected to get affected when contributions from the hyperons and 
kaons are also included in the EOS. 

In Fig. ^2 we have plotted the variation of RP~^^^ as a function of the baryon density. 
Here, R is the radius (in km) of the neutron star with IAMq and P the pressure in units 
of MeV/fm^. The solid circles denote the empirical values for the RP~^^^ at the baryon 
densities p/po = 1-0,1.5 and 2.0 taken from Ref. J6l|. In Fig. ^J we plot the ratio of 
neutron star binding energy to its gravitational mass as a function of the compactness 
parameter P = ^^, R being the radius of the neutron star with gravitational mass M. The 
solid circles represent the empirical values obtained using the relation (62| |. 

-5'bind _ 0-6/3 .„„x 

Mc"^ ~ 1-0.5/?' ^ ^ 

Finally, in Fig. El we plot the variation of moment of inertia / versus the neutron star mass. 
To obtain the values of / we have used the code written by Stergioulas J63|. The results 
presented in Fig. ^J are obtained at the Keplerian frequencies. 

V. CONCLUSIONS 

We have parameterized the GSEF containing multi density dependent terms of the form 
p*^ with z/ = 1/3, 2/3 and 1 in the local part of the Skyrme interaction. The experimental 
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data for the normal and isospin-rich nuclei are used to fit the parameters of the GSEF. 
Further, a realistic equation of state for the pure neutron matter upto high densities (~ 0.8 
fm~^) is used in the fit to ensure that the Skyrme parameters so obtained can be used to 
study the neutron star properties. For the appropriate comparison we generate a parameter 
set for the SSEF using exactly the same set of data as in the case of the GSEF. Comparing 
our results for the various quantities associated with the symmetric nuclear matter (at the 
saturation density) obtained using the parameters of GSEF and SSEF we find that the earlier 
one yields larger values for the isoscalar effective nucleon mass and for the quantity L = 3p^ 
which is directly related to the slope of the symmetry energy coefficient J. The large value 
of the isoscalar effective nucleon mass is highly desirable in order to predict appropriately 
the location of the isoscalar giant quadrupole resonance. The value of L obtained for the 
GSEF is in better agreement with the ones extracted very recently J4^ from the isospin 
diffusion data. The large value of L obtained for the parameters of the GSEF gives rise to 



M 



Ml 



the larger values for the neutron skin thickness corroborating recent predictions 

The parameter sets generated in the present work for the GSEF and SSEF are used to 
calculate the strength function for the isoscalar giant resonances and some of the key prop- 
erties for non-rotating neutron star. Calculations are performed for the strength function 
of the isoscalar giant monopole, dipole and quadrupole resonances for the ^°Zr and ^osp]-, 
nuclei. We find that the results for the isoscalar giant monopole resonance are quite similar 
for the GSEF and SSEF parameters. However, the larger isoscalar effective nucleon mass in 
the case of GSEF lowers the energy of the isoscalar giant dipole and quadrupole resonances 
by ~ 0.5 MeV and thereby improve the agreement with experimental data. The results for 
the mass-radius relationship and some of the key properties for the neutron star with the 
"canonical" mass of IAMq and for the one with the maximum mass are more or less same 
for GSEF and SSEF. We have also studied the variation of the RP~^/'^ (see Fig. El) as a 
function of the baryon density which reasonably agrees with the empirical values J61|. The 
comparison with the empirical values were also done for the variation of the ratio of the 
binding energy of the neutron star to its gravitational mass as a function of the compactness 
parameter. 



18 



Ackno^vledgments 

We would like to thank J. N. De and S. K. Samaddar for their helpful comments and for 
reading the manuscript. We also acknowledge the discussions with D. Bandyopadhyay. 



19 



[1] J. S. Bell and T. H. R. Skyrme, Phil. Mag. 1, 1055 (1956). 

[2] D. Vautherin and D. M. Brink, Pliys. Rev. C 5, 626 (1972). 

[3] L. Zamick, Phys. Lett. B45, 313 (1973). 

[4] N. V. Giai and H. Sagawa, Phys. Lett. 106B, 379 (1981). 

[5] J. Bartel, P. Quentin, M. Brack, C. Guet, and H. B. Hakansson, Nucl. Phys. A386, 79 (1982). 

[6] F. Tondeur, M. Brack, M. Farine, and J. M. Pearson, Nucl. Phys. A420, 297 (1984). 

[7] J. M. Pearson, Y. Aboussir, A. K. Dutta, R. C. Nayak, and M. F. Tondeur, Nucl. Phys. A528, 
1 (1990). 

[8] P. G. Reinhard and H. Flocard, Nucl. Phys. A584, 467 (1995). 

[9] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. Schaeffer, Nucl. Phys. A635, 231 
(1998). 

[10] B. A. Brown, Phy. Rev. C 58, 220 (1998). 

[11] M. Bender, P. H. Heenen, and P.-G. Reinhard, Rev. Mod. Phys. 75, 121 (2003). 

[12] S. Goriely, M. Samyn, P.-H. Heenen, J. M. Pearson, and M. Onsi, Nucl. Phys. A750, 425 

(2005). 

[13] B. K. Agrawal, S. Shlomo, and V. K. Au, Phys. Rev. C 72, 014310 (2005). 

[14] T. Duguet and P. Bonche, Physical Review C 67, 054308 (2003). 

[15] B. Cochet, K. Bennaceur, P. Bonche, T. Duguet, and J. Meyer, Nucl. Phys. A731, 34 (2004). 

[16] B. Cochet, K. Bennaceur, J. Meyer, P. Bonche, and T. Duguet, Int. J. Mod. Phy. E 13, 187 

(2004). 

[17] A. Bhattacharyya and R. J. Furnstahl, Nucl. Phys. A747, 268 (2005). 

[18] N. Kaiser, S. Fritsch, and W. Weise, Nucl. Phys. A697, 255 (2002). 

[19] K.-F. Liu, H. Luo, Z. Ma, and Q. Shen, Nucl. Phys. A534, 1 (1991). 

[20] K.-F. Liu, H. Luo, Z. Ma, Q. Shen, and S. A. Mozkowski, Nucl. Phys. A534, 25 (1991). 

[21] K.-F. Liu, Z. yu Ma, M. Feng, and Q. Shen, Nucl. Phys. A534, 48 (1991). 

[22] M. Farine, J. M. Pearson, and F. Tondeur, Nucl. Phys. A615, 135 (1997). 

[23] M. Beiner, H. Flocard, N. V. Giai, and P. Quentin, Nucl. Phys. A238, 29 (1975). 

[24] M. Bender, J. Dobaczewski, J. Engel, and W. Nazarewicz, Phys. Rev. C 65, 054322 (2002). 

[25] S. Shlomo, Rep. Prog. Phys. 41, 957 (1978). 

20 



[26] W. Bertozzi, J. Friar, J. Heisenberg, and J. W. Negele, Phys. Lett. B41, 408 (1972). 

[27] S. Shlomo and W. G. Love, Phys. Scr. 26, 280 (1982). 

[28] S. Goriely, M. Samyn, P.-H. Heenen, J. M. Pearson, and F. Tondeur, Phys. Rev. C 66, 024326 

(2002). 
[29] S. Kirkpatrik, J. Stat. Phys. 34, 975 (1984). 
[30] L. Ingber, Mathl. Comput. Modelling 12, 967 (1989). 
[31] B. Cohen, Master's thesis, Tel-Aviv University (1994). 
[32] T. Biirvenich, D. G. Madland, J. A. Maruhn, and P.-G. Reinhard, Phys. Rev. C 65, 044308 

(2002). 
[33] T. J. Biirvenich, D. G. Madland, and P.-G. Reinhard, Nucl. Phys. A744, 92 (2004). 
[34] G. Audi, A. H. Wapstra, and C. Thibault, Nucl. Phys. A729, 337 (2003). 
[35] E. W. Otten, in Treatise on Heavy-Ion Science, vol. 8 (ed. D. A. Bromley, Plenum, New York, 

1989). 
[36] H. D. Vries, C. W. D. Jager, and C. D. Vries, At. Data Nucl. Data Tables 36, 495 (1987). 
[37] N. Kalantar-Nayestanaki, , H. Baghaei, W. Bertozzi, S. Dixit, J. M. Finn, C. E. Hyde- Wright, 

S. Kowalski, R. W. Lourie, C. P. Sargent, et al., Phys. Rev. Lett. 60, 1707 (1988). 
[38] S. Platchkov, A. Amroun, P. Bricault, J. M. Cavedon, P. K. A. de Witt Huberts, P. Dreux, 

B. Frois, C. D. Goodman, D. Goutte, J. Martino, et al., Phys. Rev. Lett. 61, 1465 (1988). 
[39] A. Bohr and B. R. Mottelson, Nuclear Structure, vol. Vol. I (Benjamin, New York, 1969). 
[40] B. A. Brown, private communication (2005). 

[41] D. H. Youngblood, H. L. Clark, and Y. W. Lui, Phys. Rev. Lett. 82, 691 (1999). 
[42] R. Wiringa, V. Fiks, and A. Fabrocini, Phys. Rev. C 38, 1010 (1988). 
[43] P. G. Reinhard, Nucl. Phys. A649, 305c (1999). 
[44] L.-W. Chen, C. M. Ko, and B.-A. Li, Nucl-th 0509009 (2005). 
[45] M. Kutschera, Phys. Lett. B340, 1 (1994). 
[46] M. Kutschera and J. Niemiec, Phys. Rev. C 62, 025802 (2000). 

[47] J. M. Lattimer, C. J. Pethick, M.Prakash, and P. Haensel, Phys. Rev. Lett. 66, 2701 (1991). 
[48] J. R. Stone, J. C. Miller, R. Koncewicz, P. D. Stevenson, and M. R. Strayer, Phys. Rev. C 

68, 034324 (2003). 
[49] O. Bohigas, A. M. Lane, and J. Martorell, Phys. Rep. 51, 267 (1979). 
[50] B. A. Brown, Phys. Rev. Lett. 85, 5296 (2000). 

21 



[51] C. J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 5647 (2001). 

[52] A. Krasznahorkay, H. Akimune, A. van den Berg, N. Blasi, S. Brandenburg, M. Csatlos, 

M. Fujiwara, J. Gulyas, M. Harakeh, and M. H. et al, Nucl. Phys. A731, 224 (2004). 
[53] J. P. Blaizot, J. F. Berger, J. Decharge, and M. Girod, Nucl. Phys. A591, 435 (1995). 
[54] V. Kolomietz and S. Shlomo, Phys. Rev. C 61, 064302 (2000). 
[55] S. Shlomo and G. F. Bertsch, Nucl. Phys. A243, 507 (1975). 
[56] S. Shlomo and A. I. Sanzhur, Phys. Rev. C 65, 044310 (2002). 
[57] B. K. Agrawal, S. Shlomo, and V. K. Au, Phys. Rev. C 70, 014308 (2004). 
[58] T. Sil, S. Shlomo, B. K. Agrawal, and P. G. Reinhard, Phys. Rev. C (Submitted). 
[59] B. Gordon, P. Christopher, and S. Peter, Astrophysical Journal 170, 299 (1971). 
[60] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. Schaeffer, Nucl. Phys. A627, 710 

(1997). 
[61] J. M. Lattimer and B. F. Schutz, Astrophys. J. 629, 979 (2005). 
[62] J. M. Lattimer and M. Prakash, Astrophys. J. 550, 426 (2001). 
[63] N. Stergioulas and J. L. Friedman, Astrophys. J. 444, 306 (1995). 
[64] B. G. Todd-Rutel and J. Piekarewicz, Phys. Rev. Lett 95, 122501 (2005). 
[65] D. H. Youngblood, Y.-W. Lui, B. John, Y. Tokimoto, H. L. Clark, and X. Chen, Phys. Rev. 

C 69, 054312 (2004). 
[66] D. H. Youngblood, Y.-W. Lui, H. L. Clark, B. John, Y. Tokimoto, and X. Chen, Phys. Rev. 

C 69, 034315 (2004). 



22 



FIG. 1: The density dependence of the symmetry energy coefficient J(p) for the GSkI, GSkll and 
SSk forces. 

FIG. 2: The EOS of nuclear matter for the proton fractions Yp = 0, 0.25, 0.4 and 0.5 for the GSkI, 
GSkll and SSk forces. The sohd squares represent the EOS for pure neutron matter (i.e., Yp = 0) 
obtained for the UV14+UVII model \^. 

FIG. 3: Comparison of the results for the relative errors, AB = (i^^^P* — j^thyj^expt^ ^^ ^j^g values 
of binding energies for the nuclei considered in the fit to obtain the GSkI, GSkll and SSk forces. 
A=48 represents the result for ^^Ca nucleus. The result for "^^Ni nucleus is discussed in the text. 
The rms error for the total binding energy are 1.18, 1.72 and 1.28 MeV for the GSkI, GSkll and 
SSk forces, respectively. 

FIG. 4: Comparison of the results for the relative errors, Ar^h = (''c/^^ ~ ^ch)/''"ch' ' ™ ^^^ values 
for the charge rms radii for nuclei considered in the fit to obtain the GSkI, GSkll and SSk forces. 
The rms error for the charge rms radii are 0.023, 0.026 and 0.021 fm for the GSkI, GSkll and SSk 
forces, respectively. 

FIG. 5: The single particle energies for ^^^Pb obtained using GSkI, GSkll and SSk forces are 
compared with the experimental data for (a) neutrons and (b) protons. The rms error in the single 
particle energies are 1.44, 1.13 and 1.67 MeV for the GSkI, GSkll and SSk forces, respectively. 

FIG. 6: Comparison of neutron skin thikness, r„ — r^ for the GSkI, GSkll and SSk forces. The 
values for neutron skin thikness for ^^^Sn and ^ospj^ nuclei represented by filled circles are the 
recent prediction based on isospin diffusin data 44 1. A=48 represents the result for '^^Ca nucleus. 



The result for ^^Ni is discussed in the text. 

FIG. 7: Variations of the strength function with excitation energy for the (a)isoscalar giant 
monopole resonance, (b)isoscalar giant dipole resonance and (c)isoscalar giant quadrupole reso- 
nance. The results are obtained using GSkI, GSkll and SSk forces for the 208p^ nucleus. 

FIG. 8: Variations of the proton and electron fractions with the baryon density for the GSkI and 
SSk forces. 
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FIG. 9: The neutron star mass as a function of the central baryon density p^. for the GSkI and 
SSk forces. 

FIG. 10: Relation between the neutron star mass and its radius R for the GSkI and SSk forces. 

FIG. 11: Variation of RP^^''^ as a function of the baryon density p. Here, R is the radius (in 
km) for the neutron star of IAMq and P is the pressure in units of MeV/fm^. The solid circles 
represent the empirical values taken from Ref. 6ll |. 



FIG. 12: Ratio of the neutron star binding energy -Ebind to its corresponding gravitational mass 

GM 
'Re? 



M as a function of the compactness parameter (3 = %-3- . The solid circles represent the empirical 



values calculated using Eq. (|S^ . 

FIG. 13: Moment of inertia I as a function of the neutron star mass. The values of / are obtained 
at the Kepler ian frequencies. 
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TABLE I: The lower (vq) and upper (vi) limits, maximum displacement (d) and initial values 
(vj„) for the Skyrme parameters used to minimize the x^ value within the SAM. 





vo 


Vl 


d 


Vm 


to(MeV-fm3) 


-3000.0 


-1500.0 


50.0 


-1603.0 


ti(MeV-fm5) 


-500.0 


500.0 


20.0 


515.9 


t2(MeV-fm5) 


-500.0 


500.0 


20.0 


84.5 


t3i(MeV-fm3("i+i)) 


1000.0 


3000.0 


50.0 


1333.3 


t32 (MeV-fm3("2+i)) 


-1000 


0.0 


50.0 


0.0 


t33 (MeV-fm3("3+i)) 


-500.0 


500.0 


20.0 


0.0 


Xq 


-4.0 


4.0 


0.1 


-0.02 


Xi 


-4.0 


4.0 


0.1 


-0.5 


X2 


-4.0 


4.0 


0.1 


-1.713 


3^31 


-4.0 


4.0 


0.1 


0.1381 


3^32 


-4.0 


4.0 


0.1 


0.0 


3^33 


-4.0 


4.0 


0.1 


0.0 


ai 


1 

3 


1 

3 





1 

3 


a2 


2 
3 


2 
3 





2 
3 


as 


1 


1 





1 


M/o(MeV-fm5) 


100.0 


200.0 


5.0 


125.0 


V^ (MeV) 


-3.0 


0.0 


0.2 


-2.05 


V^, (MeV) 


0.0 


2.0 


0.1 


0.697 


A 


400.0 


600.0 


10.0 


485.0 


^0 


10.0 


50.0 


1.0 


28.0 
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TABLE II: Selected experimental data for the binding energies B, charge rms radii r^h, mis radii 
of valence neutron orbits r„, single-particle energies (S-P), breathing mode constrained energies Eq 
and EOS for the pure neutron matter used in the fit to determine the parameters of the Skyrme 
interaction. 



Properties 








Nuclei 






Ref. 


B 


16,24 


0, 


40,48Ca 


, 48,56,68,78^1, ^Sg^,^ 902^, 


100,132sn, 


208pb 


m 


"Teh 






160, 


40'48ca, 56Ni, 88Sr, ^^Zr, 


208pb 




[35,36] 


ry{vld^l2) 








1^0 






m 


rv{i''^h/2) 








^^Ca 






[38] 


S-P energies 








208pb 






[39,40] 


Eo 








90Zr and ^ospb 






[4JJ 


EOS 








pure neutron matter 






[42] 
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TABLE III: The values of the Skyrme parameters for GSkI, GSkll and SSk interactions obtained 
by minimizing the ^ . 





GSkI 


GSkll 


SSk 


to(MeV-fm3) 


-1855.45 


-1855.99 


-2523.52 


ti(MeV-fm5) 


397.23 


393.08 


435.00 


t2(MeV-fm5) 


264.63 


266.08 


-382.04 


t3i(MeV-fm3K+i)) 


2309.67 


2307.15 


2372.49 


t32 (MeV-fm3("2+i)) 


-449.01 


-448.28 





t33 (MeV-fm3("3+i)) 


-53.31 








Xo 


0.1180 


0.0909 


0.6835 


Xl 


-1.7586 


-0.7203 


-0.4519 


a;2 


-1.8068 


-1.8369 


-0.9214 


3^31 


0.1261 


-0.1005 


1.0508 


3^32 


-1.1881 


-0.3529 





a:^33 


-0.4594 








ai 


1 

3 


1 

3 


0.1682 


Q;2 


2 
3 


2 
3 





"3 


1 








M/o(MeV-fm5) 


169.57 


152.28 


131.98 


K, (MeV) 


-2.9944 


-2.6683 


-1.2124 


K;' (MeV) 


0.7059 


0.6191 


0.3077 


A 


538.23 


442.56 


461.43 


^0 


49.91 


47.24 


19.33 
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TABLE IV: Nuclear matter properties for the GSkI, GSkll and SSk interactions at the minimum 
value of x^- The quantities given below are: B/A the binding energy per nucleon, K^o the nuclear 
matter incompressibility coefficient, J the symmetry energy, L = 3p^ related to the slope of the 
symmetry energy, m* /m is the ratio of the isoscalar effective nucleon mass to the bare nucleon mass 
and ps the saturation density. Values for these quantities are obtained at the saturation density. 



Parameter 


GSkI 


GSkll 


SSk 


B/A (MeV) 


16.02 


16.13 


16.15 


Koo(MeV) 


230.20 


233.60 


229.17 


J (MeV) 


32.03 


34.15 


33.49 


L (MeV) 


63.46 


66.82 


52.75 


m*/m 


0.78 


0.79 


0.72 


Ps 


0.159 


0.159 


0.161 



TABLE V: The rms radii of the valence neutron orbits r^ (fm) and the breathing mode constrained 
energies (MeV). The experimental values ( and the theoretical error a) used in the fit to determine 



the Skyrme parameters are taken as follows: the values of r^ were taken from Ref. 
fm) and values of the breathing mode constrained energies were taken from |4 



3:3, 



{a = 0.06 



I ( cj = 0.5 MeV 

for the ^'^Zr nucleus and a = 0.3 MeV for ^''^Pb nucleus). The superscript 'a' denotes that the 
corresponding data were not included in the fit. 







Exp. 


GSkI 


GSkll 


SSk 


Tv 


(1^^5/2) 


3.36 


3.35 


3.37 


3.34 




(^1/7/2) 


3.99 


4.01 


4.02 


4.05 


Eo 


90Zr 


17.81 


18.12 


18.17 


18.03 




IW^^a 


15.90 


16.52 


16.59 


16.43 




144 Sm'^ 


15.25 


15.62 


15.67 


15.52 




208pb 


14.18 


13.69 


13.70 


13.66 
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TABLE VI: Values of the centroid energy ( ^ 
obtained using strength function calculated 



data given here are taken from Refs. 



41 



65, 



, constrained energy \j^^^ and scaling energy 
within the RPA frame work. The experimental 



Nuclei 


Interaction 


ma 


V ™-i 


y mi 


r 


ISGMR 


90 Zr 


GSkI 


18.48 


18.25 


19.14 


2.70 




GSkll 


18.60 


18.36 


19.25 


2.72 




SSk 


18.43 


18.20 


19.05 


2.63 




Exp. 


17.89 


17.81 






208pb 


GSkI 


13.79 


13.55 


14.46 


2.33 




GSkll 


13.79 


13.55 


14.46 


2.33 




SSk 


13.87 


13.64 


14.53 


2.31 




Exp. 


14.17 


14.18 






ISGDR 


90Zr 


GSkI 


28.56 


28.11 


29.74 


4.82 




GSkll 


28.48 


28.03 


29.69 


4.87 




SSk 


28.83 


28.38 


30.04 


4.90 




Exp. 


26.7 








208pb 


GSkI 


24.14 


23.87 


24.95 


3.60 




GSkll 


23.78 


23.50 


24.65 


3.69 




SSk 


24.61 


24.32 


25.46 


3.73 




Exp. 


22.20 








ISGQR 


90Zr 


GSkI 


14.65 


14.50 


15.52 


2.65 




GSkll 


14.42 


14.29 


15.22 


2.50 




SSk 


15.10 


14.98 


15.86 


2.50 




Exp. 


14.30 








208pb 


GSkI 


11.34 


11.23 


12.10 


2.10 




GSkll 


11.07 


10.97 


11.82 


2.05 




SSk 


11.92 


11.80 


12.71 


2.23 




Exp. 


10.89 
90 









TABLE VII: Some of the key properties of the neutron star with IAMq and the maximum mass 
Minax calculated using the GSkI and SSk forces. 







1.4M0 




-'"max 


GSM 




SSk 


GSkI 


SSk 


Pc (fm-3) 


0.53 




0.55 


1.21 


1.22 


ejc^ (1014 g cm-3) 


9.68 




10.03 


27.97 


28.54 


M{Mq) 


1.40 




1.40 


1.95 


2.00 


R (km) 


11.97 




11.62 


10.05 


9.87 


A (10^7) 


1.81 




1.82 


2.68 


2.77 


^bind (10^3 erg) 


2.03 




2.09 


5.09 


5.59 


^surf 


0.24 




0.25 


0.53 


0.58 
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